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Statement  of  Problems  Studied 

During  the  period  covered  by  this  report  we  have  studied  the  global 
dynamics  of  incompressible,  viscous  fluid  motion.  This  activity  included  a 
study  of  discrete  dynamical  systems  generated  by  temporal  and  spatial 
discretizations  of  the  Navier-Stokes  equations.  We  investigated  the 
existence  of  absorbing  sets  and  attractors  associated  with  the  two 
dimensional  Navier-Stokes  equations.  The  dynamical  properties  of  various 
finite-difference  and  finite-element  schemes  for  these  equations  were  also 
considered.  Our  goal  was  to  investigate  the  existence  of  spurious  steady- 
state  solutions.  A  resolution  of  the  issues  investigated  here  could  have 
enormous  impact  on  novel  approaches  to  the  numerical  simulation  of 
turbulent  flows. 

A  second  research  project  partially  supported  by  this  grant  included  the 
development  of  mathematical  theories,  numerical  methods,  and 
supercomputing  algorithms  for  the  control  of  fluid  flows.  Primary  areas  of 
study  (with  the  corresponding  application  of  interest)  have  been: 
Magnetohydrodynamic  (MHD)  flows  (fusion  plasma);  thermally  coupled 
mixed  fluid  flows  (chemical  vapor  deposition);  and  free  surface  flows 
(polymer  extrusion).  Results  from  the  study  of  these  problems  are 
important  to  many  Army  applications.  For  example,  this  work  could  have 
an  impact  on  the  development  of:  improved  battery  components; 
electromagnetic  armor;  and  advanced  electric  gun  propulsion  technology. 


Summary  Of  The  Most  Important  Results 

Our  initial  investigations  were  concentrated  on  the  control  of  electrically 
conducting  flows.  A  simplified  model,  involving  the  velocity,  pressure  and 
electrical  potential  fields,  was  considered.  This  model  is  valid  in  a  variety 
of  applications,  e.g.,  in  the  modelling  of  a  two-dimensional  magnetic  pump 
with  electrodes  on  the  bounding  surface,  in  the  modelling  of  liquid  lithium 
blanket  in  fusion  technology  and  metal  melt  cooling  processes.  In  the 
isothermal  case  we  have  obtained  mathematical  results  such  as  the  existence 
of  optimal  solutions  and  justification  of  the  use  of  Lagrange  multiplier 
rules  to  derive  an  optimality  system  of  equations.  We  have  also  developed 
effective  numerical  algorithms  and  successfully  implemented  these 
algorithms  in  the  case  of  simple  geometries.  See  attached  computational 
results  shown  in  Figures  MHD  a)  —  f). 
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In  addition,  we  have  obtained  results  for  control  problems  for  free  surface 
flows  in  the  context  of  extrusion  processes  for  viscoelastic  fluids.  The  two 
control  objectives  investigated  were:  (1)  to  obtain  a  quasi-uniform 
temperature  distribution;  and  (2)  to  obtain  a  desired  shape  for  the  free 
surface.  The  means  of  controls  were  heating/cooling  on  the  bounding 
surface/extrusion  die  and  adjusting  the  extrusion  speed.  Several  of  these 
results  were  reported  in  the  a  paper  on  control  in  visco-elastic  extrusion 
processes  [1],  We  also  obtained  analytical  and  computational  results  for  the 
optimal  control  of  the  steady-state  von  Karman  equations  with  distributed 
controls.  Existence  of  optimal  solutions  was  established.  Lagrange 
multipliers  were  used  to  enforce  the  constraints  and  to  derive  an  optimality 
system.  Finite  element  approximations  of  solutions  for  the  optimality 
system  were  developed  and  corresponding  error  estimates  were  derived. 
See  [2]  and  [3]. 

Finally,  we  obtained  results  on  an  optimal  control  problem  that  is  useful 
for  the  reduction  or  prevention  of  chemical  corrosion.  The  optimal  control 
problem  used  here  was  developed  for  impressed  cathodic  systems  in 
electrochemisty.  The  control  in  this  problem  was  the  current  density  on  the 
anode.  A  matching  objective  functional  was  considered.  We  first 
demonstrated  the  existence  and  uniqueness  of  solutions  for  the  governing 
partial  differential  equation  with  a  nonlinear  boundary  condition.  Next,  an 
existence  proof  of  an  optimal  solution  was  given.  Then,  we  derived  a 
necessary  condition  of  optimality.  Finally,  a  finite  element  algorithm  was 
developed  and  optimal  error  estimates  derived.  See  [4]. 
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